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Spin injection and detection in silicon is a difficult problem, in part because the weak spin-orbit 
coupling and indirect gap preclude using standard optical techniques. We propose two ways to 
overcome this difficulty, and illustrate their operation by developing a model for spin-polarized 
transport across a heterojunction. We find that equilibrium spin polarization of holes leads to a 
strong modification of the spin and charge dynamics of electrons, and we show how the symmetry 
properties of the charge current can be exploited to detect spin injection in silicon using currently 
available techniques. 



In addition to its central role in conventional electron- 
ics, silicon has spin-dependent properties (such as long 
spin relaxation and decoherence times) that could be par- 
ticularly useful in spin-based quantum-information pro- 
cessing and spintronics 0] . Unfortunately, the underlying 
origins of these attractive properties — the indirect band 
gap, weak spin-orbit coupling, and extremely small con- 
centration of paramagnetic impurities Q — also preclude 
using the standard optical methods of spin injection and 
detection in semiconductors. 

Circularly polarized light can be used to polarize carri- 
ers in semiconductors with a direct band gap. Moreover, 
both the direction and the magnitude of optically gener- 
ated charge currents 0, 0] and pure spin currents 0, Q 
can be controlled optically. In the reverse process, the 
presence of polarized carriers in a direct-gap semicon- 
ductor can be detected by measuring the circular po- 
larization of the recombination light [jj, |j| . Typical de- 
tection schemes use spin light-emitting diodes (LED's), 
where the selection rules for radiative recombination pro- 
cesses can be used to relate the circular polarization of 
the emitted light to the spin polarization of the carri- 
ers 0I1E3,E3 

For silicon, the indirect band gap makes direct appli- 
cation of these techniques problematic. We propose two 
approaches to overcome this difficulty. The first, shown 
in Fig.^a), is based on high-quality heterojunctions be- 
tween Si and a direct-gap semiconductor. In such a het- 
erojunction, optical techniques could be readily employed 
in the direct-gap semiconductor to circumvent the prob- 
lems with spin injection and detection in Si. The key pre- 
requisite for such a proposal is an interface with Si that 
would not be detrimental to the spin transport. This is a 
nontrivial undertaking, as the lattice mismatch between 
Si and most direct-gap semiconductors typically leads to 
low-quality interfaces with a high density of interfacial 
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FIG. 1: Proposed schemes for spin injection and detection in 
silicon, (a) Optical realization based on radiative processes 
(excitation for spin injection and recombination for detection) 
in direct-gap semiconductors surrounding silicon. Arrows de- 
pict the spatial decay of nonequilibrium spin, (b) Electrical 
realization based on spin splitting and net spin density (mag- 
netization) in the two magnetic regions. The relative orienta- 
tion of the nonequilibrium spin in Si and the equilibrium spin 
in the magnetic regions influences the magnitude of a charge 
current or an open-circuit voltage. Other realizations are also 
possible by combining schemes (a) and (b). 



defects. Nevertheless, there has been recent progress in 
fabricating high-quality GaAs/Si interfaces (despite the 
4% lattice mismatch) |l2j, and GaPi-^N^ leads to even 
smaller mismatch (below 1%) Recent studies of 

charge transport in GaAs/Si heterojunctions suggest the 
feasibility of the scheme shown in Fig.^a). In particular, 
GaAs/Si heterojunctions displayed I-V characteristics of 
an ideal diode |14j . and optical excitations were studied 
in GaAs LED's grown on Si 

In the second approach, shown in Fig. ^b), mag- 
netic semiconductors approximately lattice matched with 
Si could be used for spin injection and detection [lfij . 
For example, the Mn-doped chalcopyrite ZnGeP2 (mis- 
match < 2%) mm 

has been reported to be ferromag- 
netic at room temperature. Another Mn-doped chalcopy- 
rite, ZnSiP2, was recently predicted m to be ferromag- 
netic, as well as highly spin polarized and closely lattice- 
matched with Si (mismatch < 1%). Mn doping of the 
chalcopyrite alloy ZnGei_ :E Si x P2 would likely lead to an 
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exact lattice match, since the lattice constant of Si is 
between those of closely matched ZnSiP2 and ZnGeP2- 

To capture the main features of spin-polarized trans- 
port across a hetero junction we formulate here a model, 
representing both schemes in Fig. ^ that can be solved 
analytically. One can model the right (injecting) elec- 
trode by appropriate boundary conditions, and hence we 
focus on the two left regions that define the heterojunc- 
tion in our model. A heterojunction is doped with ionized 
acceptors and donors of density N a and N4. Codoping 
with magnetic impurities would additionally introduce a 
net spin, but need not change the number of carriers. An 
inhomogeneous distribution of N at( i implies a large devia- 
tion from local charge neutrality, so that Poisson's equa- 
tion must be explicitly solved. For nondegenerate doping 
levels (Boltzmann statistics) the spin-resolved quasiequi- 
librium electron and hole densities are 

(1) 

where A = +1 for spin up (|) and —1 for spin down (J.). 
The total electron density n = n-\ + ni can also be de- 
composed as a sum of equilibrium and noncquilibrium 
parts, n = no + 5n. We define the electron spin den- 
sity s n = n-f — 7i| and the spin polarization P n — s n /n, 
with an analogous notation for holes. In Eq. QJ, sub- 
scripts c and v label quantities pertaining to conduction 
and valence bands. The corresponding effective density 
of states are N CtV — 2{2itm c ^ v kBT /h 2 ) 3 / 2 , where m c> „ are 
effective masses. The spin-A conduction band edge (see 
Fig. |5J) E c \ = E c0 — q<j) — Xq( c differs from its nonmag- 
netic bulk value E c q because of the electrostatic poten- 
tial <j) and the spin splitting 2q(, c , which parameterizes 
Zeeman or exchange splitting due to magnetic impurities 
and/or an applied magnetic field 0. Here, AE C is the 
conduction band edge discontinuity and [i n \ — /i + 5\i n \ 
is the chemical potential for spin-A electrons. An analo- 
gous notation holds for the valence band and holes. 

We assume transport across the interface is domi- 
nated by drift-diffusion, so that the spin-resolved charge- 
current densities are 

J«a = qfi n xn x \7E cX + qD nX N c V(n x /N c ), (2) 
J P A = qjip\P\VE v x - qDp X N v \7(p\/N v ), (3) 

where /2 and D are mobility and diffusion coefficients. 
We note that "drift terms" have quasi-electric fields 
<x VE C>V \ that are generally spin-dependent (VCc,u 7^ 
is referred to as a magnetic drift |20|) and different for 
conduction and valence bands. In contrast to homojunc- 
tions, additional "diffusive terms" arise due to the spatial 
dependence of m CjV , and therefore of N c>v . 
We write the continuity equation for n\ as 

-dnx/dt + V ■ J n \/q = +rx{n x p\ - nxopxo) (4) 
+ [n x - n_ A - XP n on]/2T sn - G\, 
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FIG. 2: Band diagram for a magnetic heterojunction. In equi- 
librium, the chemical potential no is constant. Conduction- 
and valence-band edges (E c and E v ) are spin-split in the mag- 
netic p-region, while there is no spin splitting in the nonmag- 
netic n-region (corresponding to Si). For a sharp doping pro- 
file, there are generally discontinuities (AE C and AE V ) in the 
conduction and valence bands at x — w/2. 



with an analogous equation for p\. Here, r\ is the recom- 
bination rate of spin-A carriers; T sn . P is spin relaxation 
time for electrons and holes; and G\ is the photoexcita- 
tion rate due to electron-hole pair generation and optical 
orientation (when Gj ^ Gj.) 0-0] Spin relaxation equi- 
librates carrier spin while preserving noncquilibrium car- 
rier density , so that for nondegenerate semiconductors 
we have P n o — tanh(<7£ c /fcBT). 

We make several assumptions allowing us to solve this 
model analytically. We focus here on the steady-state 
low-injection regime, at applied bias \V\ < min{\E c \ — 
E v \\). Spin-orbit coupling in the valence band typically 
leads to a much faster spin relaxation of holes than elec- 
trons (3-4 orders faster in GaAs Q), and so it is reason- 
able to consider that the spin of holes is in equilibrium 
(Ss p — and P p — P p o). As a result, only the noncqui- 
librium electron spin density (5s n — ► 5s) and the mi- 
nority carrier density need to be calculated throughout 
the heterojunction. We assume a sharp doping profile 
Nd(x) — N a (x); referring to Fig.|2| this leads to a discon- 
tinuous change in materials parameters at x — w/2. We 
take N c>v , p,, D, and the permittivity e to be constant 
outside the space-charge region xl < x < xr, and hence 
label them by indices L and R. The width of a space- 
charge region is xr~xl oc (Vm — V) 1 ^ 2 , where the built-in 
voltage is qV bl = -AE C + k B T\n(nQ R N cR /n 0L N cL )\ we 
note that the discontinuity AE CjV can be accurately mea- 
sured at interfaces with Si [ilj . Equations @ and © , to- 
gether with the continuity equations, reduce to diffusion- 
like equations for Sn, 5s in the p-region and 5p, 5s in the 
n-region. For the the (magnetic) p-region, we find that 
the spatial dependence of both 5n and 5s are described 
by two distinct decay lengths; this is in marked contrast 
to previously studied cases These decay lengths can 
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FIG. 3: Decay lengths k _1 and x _1 in the magnetic p-region, 
normalized to the electron diffusion length L n = (Z?nT„) 1//2 , 
as a function of hole-spin polarization for fixed electron-spin 
polarization P n o — 0.35. Each curve is labeled by the ratio of 
electron lifetime to spin-relaxation time, r n /r s „. 



be written as 



[(L- 2 +LJ 2 + L-*)/2 



(5) 
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p n0 L; 2 ) 2 + (i-p„ 2 )i; 4 ] 1/2 /2]- 1/2 



where the upper sign refers to n~ l ; L\ = [p n /r\p\a) 1 / 2 
are the electron diffusion lengths; and L s = (DnTsn) 1 / 2 
is the electron spin diffusion length. 

It is instructive to consider the regime of spin- 
unpolarizcd holes, appropriate for the scheme shown in 
Fig. QJa). In this regime we have pj = Pio = N a /2 
and r-f = rj, = r/2 (in a nondegenerate regime r-f « r*i 
even for P p ^ [13 )• It follows from Eq. © that k _1 
reduces to L sn — [D n T a ) x / 2 , where i sn is the effective 
electron-spin diffusion length and T s = (rN a + l/r^)" 1 
is the electron-spin lifetime. Analogously, x" 1 reduces 
to the electron diffusion length L n = JDnTn) 1 / 2 , where 
r n = l/rN a is the electron lifetime Thus, in this 
regime k _1 and x _1 separately determine the decay 
lengths for Ss and Sn, respectively. For the more gen- 
eral case, in Fig. [21 we show k^ 1 and x -1 as function of 
P p . From the behavior of the decay lengths for Sn and 
Ss, we conclude that polarization of the equilibrium hole 
spins leads to a strong modification of charge and spin 
dynamics of electrons. 

We turn now to the more general case of spin detec- 
tion using magnetic semiconductors, shown in Fig.Q](b), 
and solve the corresponding problem of spin-polarized 
transport across the heterojunction in Fig. [21 We im- 
pose the ohmic boundary conditions Sn = Ss = at 
x = 0, and include optical or electrical carrier and spin 
injection through the boundary conditions Sp ^ and 
Ss 7^ at x = w. To match the chemical potentials 
Hn,p\ at xl and xr (which is an accurate approxima- 
tion of a full numerical solution in magnetic p-n junc- 
tions |2ojp requires satisfying the self-consistency condi- 



tion = [P& + SP*)/(l + P£ SP?), where SP* = 
Ssn/Nd is determined from the continuity of the spin cur- 
rent, D n Ld(SsL)/dx — D n jid(5sn)/dx. Consistent with 
this matching, a generalization of Shockley's relation j2^| 
is 5n L = n 0L [e^p(qV/k B T) - 1] + s QL exp(qV/k B T)SP^ 
and(5s L = s 0L [exp(qV/k B T)-l]+noLeiq>(qV/kBT)SpK, 
such that 5nL,SsL and 5nR,SsR can be considered as 
boundary conditions in the p- and n-region, respectively. 
For the p-region we then obtain 



Sn, Ss = 



where C K 



(X 2 + a)Sn L + b5s L 
(k 2 — x 2 ) sinh(«;a;L) 
(k 2 + a)SnL + bSsL 
' (k 2 — x 2 ) sinh(xcci) 



sinh(«a;) (6) 



sinh(xs), 



1 for Sn; C KS 



-(k 2 + a)/b and 



C xs — — (x 2 + a)/b for Ss; and we have defined a = 
-(£f 2 + LJ 2 )/2 and b = ~{L^ 2 - LJ 2 )/2. 

In the nonmagnetic n-region (representing Si), the to- 
tal charge current J is the sum of minority carrier cur- 
rents at xl and xr^ J — J u l + J p r, in analogy to Shock- 
ley's formulation |24j. Here, J n L — qD n i,d(Sni,)/dx, and 
d(Sn L )/dx can be evaluated using Eq. ©. A straight- 
forward method for detecting injected spin in Si follows 
from the symmetry properties of the different contribu- 
tions to the charge current under magnetization reversal. 
By reversing the equilibrium spin polarization using a 
modest external magnetic field (P n o,P p — > —P n a,—P p ) 
it follows from Eq. I© that k, x — ► X an d a,b —* a, —b. 
A part of JnL, odd under such reversal, can be identified 
as the spin- voltaic current |2o| . 

J sv cx qD nL n 0L e^{qV/k B T)SP^/{n 2 - X 2 ), (7) 

which originates from the interplay of the equilib- 
rium and nonequilibrium (injected) spin polarization in 
the p- and the n-region, respectively. Measurements 
of J(V,P n0 ,P p ) - J(V,-P n0 ,-Pp) = 2J sv (V,P nQ ,P P ) 
would then provide: (1) cancellation of contributions to 
the charge current that are not related to the injected 
spin in Si; (2) a choice of V to facilitate a sufficiently 
large J sv for accurate detection. Alternatively, for the 
optical detection scheme of Fig. ^a), one would consider 
the limit P n0 = P p = in Eq. ijfj}, and evaluate P n {x). 

We illustrate the effects of spin injection and detection 
in Si using the heterojunction in Fig. [3 and show the 
results in Fig. 0] We use a standard set of parameters for 
Si doped with N d = 10 17 cm" 3 25]: N cR = 3.2 x 10 19 
cm" 3 , D nR = AD pR = 20 cmVa, e R = 11.7, t p = 10" 7 s, 
an intrinsic carrier density |23j mn = 10 10 cm -3 , and es- 
timated spin relaxation time in Si as t s r = 10~ 7 s. Circu- 
larly polarized light in zinc-blende semiconductors, such 
as GaAs, generates spin polarization of up to SP n = 0.5, 
which can be increased even further using strain or quan- 
tum confinement 0, 0] • We model the effects of spin in- 
jection from a neighboring direct-gap semiconductor, as 
depicted in Fig. ^ by assuming 5P n = 0.4 at x = w, and 
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FIG. 4: Electron-spin polarization, P n , for a magnetic hetero- 
junction at forward bias V=0.7 V. The equilibrium spin po- 
larizations in the p-region are P„o =0.35 and P p = 0.8, while 
the injected spin polarization in the n-region is SP n = 0.4 at 
w — 3 /jm. Inset (a): profile of the electron carrier and spin 
density in the p-region. Inset (b): voltage dependence of the 
electron charge current J n and its part due to the nonequilib- 
rium spin (the spin- voltaic current J av ). 
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10, r n = 10- y s, n lL = 10 b 
cm - ", and r s ^ = 10~ 10 s. We parameterize the spin 
splitting of carrier bands (see Fig.|2J) with P n0 — 0.35 and 
P p = 0.8 which, for Boltzmann statistics, corresponds to 
a splitting of q( v w 3q( c ss /c^T [5(J. The effect of decay 
lengths in the p-region (k^ 1 rs 0.22 /zm and x _1 w 0.64 
/im) can be seen in the Fig. Ufa). In the n-region, 
is/} = (DnBjsii) 1 / 2 « 14 pm implies a weak spin decay 
(6P™ ~ 8Pn)i an d the self-consistency condition gives an 
estimate P% w (0.35 + 0.4)/(l + 0.35 x 0.4) w 0.66. We 
find a large "signal-to-background ratio," J sv / J ni over a 
range of bias voltages [see Fig.^Jb)], and thus we predict 
that the magnitude of this spin- voltaic current — a finger- 
print for spin injection and detection in Si — could readily 
be detected using existing experimental techniques |27|. 

Our theoretical framework could be also used to calcu- 
late the effects of nonequilibrium spin in other magnetic 
hetero junctions, including Si-based spin transistors |28| . 
An extension of our results to the regime of large for- 
ward bias could provide an alternative method to de- 
tect nonequilibrium spin in Si. Even without a source 
of spin injection, a nonequilibrium spin could be gen- 
erated through the process of spin extraction [2(1 l2fj| 
into a neighboring magnetic region. Our findings should 
also provide useful test cases for developing more gen- 
eral numerical methods for treating spin-polarized trans- 
port [il. 
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